Spin-Seebeck effect in a strongly interacting Fermi gas 
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We study the spin-Seebeck effect in a strongly interacting, two-component Fermi gas and propose 
an experiment to measure this effect by relatively displacing spin up and spin down atomic clouds 
in a trap using spin-dependent temperature gradients. We compute the spin-Seebeck coefficient and 
related spin-heat transport coefficients as functions of temperature and interaction strength. We 
find that when the inter-spin scattering length becomes larger than the Fermi wavelength, the spin- 
Seebeck coefficient changes sign as a function of temperature, and hence so does the direction of the 
spin-separation. We compute this zero-crossing temperature as a function of interaction strength 
and in particular in the unitary limit for the inter-spin scattering. 



Introduction- Spin caloritronics, the study of coupled 
spin and heat transport, is a rapidly developing subfield 
of spintronics [T]. In particular, the spin-dependent gen- 
eralization of the Seebeck effect, called the spin-Seebeck 
effect, has been intensively studied in the solid-state en- 
vironment P^. Recently, there has been broad interest 
in exploring spintronic phenomena in cold atomic sys- 
tems [3H5]. Spin transport in a strongly interacting, two 
component Fermi gas was investigated experimentally in 
Ref. [6 . It is the purpose of this Letter to study the 
associated heat transport, i.e., thermo-spin effects, in a 
similar setting. 

In the ordinary Seebeck effect in metals, an electro- 
chemical potential gradient is generated by applying a 
temperature gradient. Similarly, for a gas with two spin 
states, the spin-Seebeck coefficient Ss determines the spin 
chemical potential /x^ generated by a spin temperature 
gradient VTg through the relation V/i^ = SgS/Ts, where 
fis = ^^ — /i_ and Ts = T+ —T^, fi^ and T^ being the 
spin-dependent chemical potential and temperatures of 
the spin a atoms, respectively, and we label spin compo- 
nents by -I- and — . To measure the spin-Seebeck coeffi- 
cient, we propose relatively displacing the center of mass 
of spin up and spin down atom clouds in a harmonic trap 
by applying a spin-dependent temperature gradient, for 
example by selectively heating one spin component with a 
laser, as illustrated in Fig.fl] The locations x± of the cen- 
ter of mass of the spin up and down atoms are shifted to 
the minimum of /^± -(- V ^ where V trapping potential, re- 
sulting in a spin separation Xs — x+ — X- — Ss'S/Ts/muj'^ , 
where m is the mass of the atoms, w is the trap frequency 
in the direction of the temperature gradients. For an or- 
der of magnitude estimate, we take Sg — .Olfcs, as ver- 
ified below. For VT^ = IQ-^K/cm [7], w = 27r x 1.46, 
we find x^ — 1 mm, which is well within experimental 
resolution. 

We have computed the spin-Seebeck coefficient for a 
two-component Fermi gas, plotted in Fig. [2] as a function 
of temperature and for several values of the interaction 
strength kpa, where kp is the Fermi wave vector and 
a is the inter-spin scattering length. As seen from the 
figure, for weak interactions {kpo- < 1), Ss is small and 
negative, while for strong interactions {kpa > 1), Sg is 
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FIG. 1: Spin up and spin down atomic clouds are spatially 
relatively displaced in the presence of a spin temperature gra- 
dient. The distance between the center of mass of the differ- 
ent spin components, denoted by Xs, is proportional to the 
spin-Seebeck coefficient. 



larger and its sign changes as a function of temperature. 
In terms of the experiment mentioned above, this means 
that the spin displacement changes direction as a func- 
tion of temperature, which is an interesting qualitative 
effect. We also plot the zeros of Sg as a function of kpa 
in Fig. [SH. The temperature of the zero-crossing reaches 
a universal value Tg, in the unitary limit for the inter- 
spin scattering length kpa — )■ oo. We find Tq ~ .378 Tp 
in our calculation, where Tp is the Fermi temperature. 

The thermodynamic reciprocal of the spin-Seebeck ef- 
fect is the spin-Peltier effect, in which a spin-dependent 
heat current proportional to the spin-Seebeck coefficient 
is induced by a spin current. This effect will heat up spin- 
up and spin-down components differently and provides 
another way to measure Ss- Furthermore, as discussed 
in Ref. [5], the spin-Seebeck effect contributes to the to- 
tal dissipation so that Ss can also be measured through 
the heating. 

We note that the spin-Seebeck coefficient was calcu- 
lated for a weakly interacting Bose gas in Ref. "F , but the 
Bose gas is unstable towards the formation of molecules 
for large scattering lengths, which makes the strongly in- 
teracting regime more difficult to realize experimentally. 
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FIG. 2: The Seebeck coefficient plotted as a function of the 
reduced temperature T/Tf, where Tp is the Fermi temper- 
ature, for values of kpa representing weak {kpa — .1) and 
strong (kpa = 1) coupling, and approaching the unitary limit 
{kpa = 10). 



Phenomenology- We are specifically interested in phe- 
nomena due to spindrag, the transfer of momentum be- 
tween different spins due to inter-spin scattering, which 
allows one to generated currents in one spin species by 
applying forces on the other, and we define a set of spin- 
heat transport coefficients which captures these effects as 
follows. We consider a Fermi gas with two different spin 
states selected from a larger half-integer spin multiplet, 
which we will call "spin up" (-I-) and "spin down" (— ), 
for equal spin up/down densities n^ — n_ = n, i.e., in 
the absence of spin polarization, and apply equal and op- 
posite forces and temperature gradients for the two spin 



species. I.e., 



-F_ and VT+ 



-Vr_. In linear 



response, the ensuing spin current and spin heat current 
defined by js = j+ - j-, and Qs = q+ - q-, respectively, 
are given by 



1 Ss 



F, 
-VT, 



(1) 



where F^ = F+ - F_ is the spin force, T^ ^ T+ - T_ 
is the spin temperature, T is the equilibrium temper- 
ature, CTs is the spin conductivity. Kg is the spin heat 
conductivity (at zero spin current), ZgT ~ asS^T/Ks, 
and Onsager reciprocity is explicitly included in the ma- 
trix above. We note that F^ is the thermodynamic force 
which includes forces coming from pressure gradients, 
i.e., Fs = fg'^* — 'Vps/n, where i°^^ is the external spin 
force, and Ps — p+ — P- is the difference in pressures of 
the spin up and down atoms. These coefficients, com- 
puted with the Boltzmann equation described below, are 
plotted in Fig. Is] as functions of T jTp for several values 
of kpa. 

As is well known, the spin conductivity as rapidly in- 
creases at low temperatures due to Pauli blocking. Our 
result for Og., plotted in Fig. [Sk includes corrections due 
to spin-heat coupling, but they are negligibly small, so 
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FIG. 3: (a) A plot of the spin conductivity normalized 
as fiA(Ts, (b) the spin heat conductivity normalized as 
hA.His/kB^P; and (c) the figure of merit Z^T, for kpa — .1 
(thick blue), kpa — 1 (dashed purple), kpa = 10 (dotted yel- 
low) . (d) A plot of the temperatures To relative to Tp where 
Ss = as a function of kpa. 



that one can safely take CTs — nTs/m with Ts the spindrag 
relaxation time [5| measured in Ref. j6j and calculated in 
Ref. fSl. The downturn of Sg at low temperatures is a 
quantum mechanical effect that also occurs for bosons, 
where in contrast to fermions, the spin conductivity de- 
creases sharply at low temperatures due to bosonic en- 
hancement of scattering |10j . The spin heat conductiv- 
ity Kg is plotted in Fig. [Sb, where it is seen to increase 
with increasing T. The dimensionless figure of merit ZgT 
(Fig. Isb) determines the thermodynamic efficiency of en- 
gines based on thermo-spin effects [llj . 

A spin-dependent temperature gradient can only be es- 
tablished when intra-spin scattering is much greater than 
inter-spin scattering. For fermions, the first nonvanish- 
ing intra-spin scattering amplitude is p-wave. To make 
this large, one can tune the p-wave scattering length by 
a Feshbach resonance [H]. Taking the unitary limit for 
intra-spin scattering, the intra and inter-spin differential 
cross section are given by 
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where Pr — |pr| is the relative momentum of incoming 
particles with momenta pi and p2, defined by Pr = (pi — 
P2)/2, the hat superscripts denotes unit vectors, and a 
is the inter-spin s-wave scattering length. 

Calculation of transport coefficients.- Next, we present 
the computation of Ss using the Boltzmann equation. 



We parametrize the non-equilibrium, steady-state distri- 
bution by 



npa{r) = fpair) - dj°<l)p^{r), 



(3) 



where /p = (cxp[(ep — ii)/kBT] + 1) -"^ is the equiUb- 
rium Fermi distribution, /i is the chemical potential, ep = 
p2/2m, and /p^(r,t) = (exp [tp - ^^(r))/fci3T^(r)] + 
1)^^ is the local equilibrium distribution, d^f^ = 
~/p(l ~ fp)/^BT, and 4'pcr is determined by solving 
the Boltzmann equation for the spin distribution Ups = 
Up^ — ?ip_ in linear response. 
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(4) 



where w{T) = /i + Ts is the enthalpy per particle and s 
is the entropy per particle [13j . We defined (t>ps = fc_B0p • 
(— VTs), and expressed the linearized collision integral 
in the Boltzmann equation for the spin distribution as 
{dnps/dt)co\\ = Cp[</>p] • (— VTs). The spin current is 
given by 
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(5) 



We solve Eq. Q using the method described in Ref. [5] . 
Applying the temperature gradient along the x-axis, we 
parametrize the response by a power series, 



(t)ps{a,T) = 



6o(a,T) + 6i(a,T)(^ 



Pxi-kBd^Ts) . 

(6) 

The coefRcients bo and 6i, determined by the approxi- 
mate solution to Eq. Q, are given by 



bo{a,T) 
bi{a,T) 



3nl{T) 



where 
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(8) 



z = e''/'=«^ is the fugacity, /„(z) = -Li„(-z), Li„(z) 
are the polylogarithmic functions, and C„„i are the 2x2 
matrix elements of the collision integral. 
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FIG. 4: Plots of the transport coefficients in the high- 
temperature limit, T ^ Tf, as a function of a/A. (a) The 
spin conductivity normafized as hAas, (b) the spin heat con- 
ductivity normalized as mls!^Ks/2-RhkB, (c) Ss/ks and (d) 
ZsT. 



where we define A_| [4>p\ ~ ^p, -I- (j)p^ — (pp^ — (jjp^ and 

A_| [(j)p\ = (/)p3 — (J3p^ — (pp-^ + 0p2 for an arbitrary func- 
tion 0p, and in the integrand momentum conservation is 
satisfied: Pi + P2 = P3 + P4- 

From Eq. ^ and Eq. (fTl), it follows that the See- 
beck coefficient is given in terms of fcoj^i and the spin 
conductivity (Tg by 



Ss{a,T) 



bo{a,T) + bi{a,T) 



knT 



(10) 



The other transport coefficients are calculated similarly. 
In order to make the numerics more tractable, we have 
omitted the angular dependence in dcr_(--|_/df7. On the 
other hand, in the high-temperature limit, the integrals 
in Eq. ^ including the angular factor can be done 
analytically [13], which allows us to calculate the co- 
efficients in the high-temperature limit. The result is 
shown in Fig. H plotted as a function of a/ A, where 
A = \J 2tt h^ /rakBT. In this limit, the temperature de- 
pendence of (Ts can be understood from classical consid- 
erations. The spin conductivity is approximately related 

to the inter-spin collision time t_| by as oc ht-i and 

1/t-^ = na^ vt, where a^ is the inter-spin cross sec- 
tion and vj- ex 1/A is the average thermal velocity. In 

the limit a/A -^ 0, a^ oc a?, thus Acts cx (A/a)^, and 

when a/ A — > 00, (t_| oc A^ thus Act^ oc 1, in agreement 

with our result. 

The behavior of the transport coefficients depends 
crucially on the shape of the perturbed spin distribu- 
tion 6nps — —defp<pps and the associated spin current 
density, which we plot for kpa = 10 in Fig. [5) The 
positive (negative) parts of (Jrip^. may be regarded as 
particle (holes) having group velocities ±p/m. Since 
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FIG. 5: A plot of (a) the perturbed spin distribution per spin 
temperature gradient, normalized as Sus / {—dxT / kpT) along 
the paj-axis, and a plot of (b) the spin current density along 
the j)a;-axis, per de per spin temperature gradient, normalized 
as {djs/de){h/kp(dxTs/T)), for kpa — 10 and temperatures at 
which Sa is negative (T = .3Tf, blue thick line) and positive 
(T = .5Tf, purple dashed line). 



Sripc^ = — (5np_o., every spin up particle is matched with 
a spin down hole with the same momentum, resulting 
in the spin current. Thus, the sign of Ss is determined 
by the relative number of particles or holes induced in 
response to the spin temperature gradient. 

Next, we give a criterion that determines the sign of the 
spin-Seebeck coefficient and show that it does not depend 
on the form of the intra-spin scattering at all. First, we 
note that in our solution given by Eq. ([7|, we always 
have 6o < and 6i > because l(T) > 0, Cnm > 0, and 
det C > [fSj, which implies 6o < and 5i > 0. Thus for 
positive momenta, 6i(6o) corresponds to particles (holes) 
created above (below) the Fermi surface. Inspecting Eq. 
( 10 1, we find that the criteria to have 5*8 < is 
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and the opposite inequality for Ss > 0. Thus Ss is pos- 
itive when bi becomes large enough to violate Eq. (11) 
|16) . Furthermore, it turns out that the only integral 
in Eq. ([9]) containing intra-spin scattering that is non- 
vanishing is Cii, which does not enter in Eq. ( 11 ). 

The temperature dependence of bo, bi follows from the 
temperature dependence of the collision matrix elements 
in Eq. ([9]), which are given by integrals nonvanishing only 
for pr ~ ^/4nh/A. Therefore, it is useful to express the 
differential cross section Eq. (pi) in terms of the rescaled 
momentum p^ = {A/\/4nh)pr, 
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From Eq. ( 12 ) we see that the change in the sign of Ss 



is related to the crossover from hard-sphere scattering, 
da/dfl ~ a^ when kpa <C 1 or T/Tp ^ 1, to momentum- 
dependent scattering, da/d^ ^ pr^ when kpa ~ 1 and 
T~Tp. 

Discussions and outlook- We note that as the temper- 
ature is lowered, one expects to enter the Fermi-liquid 
regime, for Tc < T <C Tf, where Tc is temperature for 
the superfluid transition, and one should use Fermi-liquid 
scattering amplitudes in the collision integral [8]. Near 
and above Tc, one should also take into account effects 
of pairing correlations on the inter-spin interaction |17) . 
Both these regimes can be analyzed with the formalism 
presented in this paper, and will be relegated to future 
work. 
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